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The chemical ordering transition in a binary alloy is examined using classical density functional theory 
for a binary mixture. The ordered lattice is assumed to be obtained from the disordered lattice by 
a volume change only as in LI2 ordering from an face centered cubic chemically disordered crystal. 
Using the simplest possible approach, second order truncation of the expansion, non-overlapping Gaussian 
distributions at the sites, and expansion of the correlation functions about the sites, a very tractable 
expansion is obtained. Under these assumptions the expansion consists of the same terms as the lattice 
gas formalism where the lattice is implicitly taken as fixed, plus additional interaction terms, and an 
additional entropy term. This additional entropy term represents a lowest order approximation to the 
vibrational entropy change. 
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I. INTRODUCTION. 



Chemical ordering transitions in alloys have been studied by a variety of methods, including a lattice gas analogue 
of classical density functional theory . Since lattice changes also occur on ordering it would be valuable to be able 
to include them along with the chemical changes in a similar approach. One approach would involve a set of variables 
including a global elastic strain tensor along with chemical occupation variables H . A question that arises in such an 
approach is the appropriate form for the ideal (non-interacting) free energy. It has been shown that the occupation 
variables alone can be treated as a complete system, and that the ideal free energy is the ideal mixing entropy. The 
occupation variables, plus a global strain tensor, are not a complete system, however. One consistent set of variables 
is the position and momentum variables for a binary mixture. The approach taken here to developing a lattice gas 
plus strain tensor formulation of ordering is to look to the classical density functional theory of the binary mixture, 
and simplify it to the lattice. Since the part of the strain tensor that most directly affects the entropy is the volume, 
the other terms will be ignored here. 

Several forms of classical density functional theory have been used to study the freezing of binary liquids. Since 
the concern here is the "entropy" terms, the simplest form of the theory will be used. Here instead of looking at the 
freezing of a liquid, the free energies of two solid structures (the disordered, and the (partially) ordered) are compared. 
A similar approach has been used by Sengupta, Krishnamurthy and Ramakrishnan to study the fcc-bcc interface || , 
and a form of density functional theory has been used to study the ordering of hard sphere mixtures M . 



II. CLASSICAL DFT FOR MULTIPLE SPECIES. 



Classical density functional theory of mixtures has been used to study the freezing of binary hard sphere and 
Lennard- Jones fluids fq-llq]. 

' Consider a system with a fixed volume V, and m species of classical particles. [The primary interest herein is 
^ , to = 2.] The partition function is computed in the grand canonical ensemble at fixed values of V, the temperature, T, 
J^j ' and the chemical potentials, [i a , a = 1, . . . , m. 

Classical DFT guarantees the existence of a functional f2 of the average densities p a (r) which is minimized by the 
■ equilibrium average densities, and which evaluated at the equilibrium densities is equal to the grand potential. For a 
non-interacting system this can be computed explicitly and is B: 

tt ld eal[{p a }] = /T 1 V / df Pa (f) (\n(X 3 aPa (f)) - /3p a - 1) , (1) 
a J V 

where j3 = l/{ksT), and A Q is the de Broglie thermal wavelength of species a. Note that equation ([!]) is simply the 
sum over the species of the ideal free energy of each species. 
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Let V$[{/9 Q }] = Clideai [{Pa}] and expand $ about some particular uniform (liquid) state, with densities 

{Pao}- Truncating at second order in the difference in densities, Sp a (see ^] |l3| for more careful alternatives), and 
letting N a = J v dr p a (f) one obtains as the expansion of the mixed functional for the grand potential: 

pn[{ Pa }} = v<f[{ Pa0 }} 

+ W drp a (r) (In (X 3 a p a (r)) - l) 

- y^^/3p a N a 



+ \ X] / df i df 2 C^(ri, r 2 ) (p a (rl) - p a0 ) (ppfa) - ftao) 

+ 0(<5p 3 J. (2) 
Multiplying out the products, and collecting constant and linear terms, this can be rewritten as: 

a Jv 
+ Vf 
+ D a N a 

a 

+ J dridr 2 C^(r l ,r 2 )p a (ri)pf}(r2), (3) 

a (3 ^ 

where 

/o = $ [K oll-E^Wo + iE^o/ dr 2 (0, r 2 ) (4) 
A* - -/3 Pa + M A 3 ) - 1 + C« - ^ P/3 o / dr 2 <7<g (0, f 2 ) (5) 

/3 

By assuming that the uniform fluid at a particular total density and set of concentrations is a local minimum of 17, 
D a could be eliminated in favor of the liquid densities {p a o}- However, here the task is to compare trial ordered 
solids with the disordered solid at temperatures well below the melting point. The D a will therefore by retained for 
now, and later eliminated using the disordered solid as the reference state. 

The constant /o drops out of the difference in grand potential between two states, and so is irrelevant to determining 
the transition state and temperature. It can be determined in terms of the pressure of the reference disordered state. 

C^ (fi, r 2 ) is symmetric in a/3, and depends only on r^— r\. (Note that for fixed {p a } these liquid direct-correlation 
functions are needed as a function of temperature in a range including the ordering temperature. In many systems 
this would involve extrapolation beyond reasonable temperatures for the liquid, and so this computation would not 
be feasible as written. The purpose here, however, puts more emphasis on understanding the appropriate form of the 
expansion for the solid, than on a computation using actual liquid data. Given the appropriate form of the expansion 
stated in terms of the occupation variables and strain tensor, other sources for the parameters are likely to be more 
suitable to actual numerical work.) 



III. DENSITY ANSATZ FOR DISORDERED AND ORDERED SOLID 

Consider a fixed volume V. A solid is described in terms of a set of real-space lattice vectors {Ri}, with N sites in 
the volume V. Making the simplifying assumption that there is exactly one atom per site, N = ^ Q N a . In general £1 

depends on {Ri}, and an assumption needs to be made as to what possible lattices are being considered. To keep the 
situation as simple as possible, assume that both the disordered solid and the trial states for the ordered solid have 



2 



the same lattice structure, so that the {Ri} depend only on N. (Up to an arbitrary choice of origin for the lattice.) 
So either TV, or the lattice constant, or the average total density can be considered as one parameter, and there is 
one constraint on the vector {N a }. Think of p = N/V as the parameter, but N will be written for pV whenever 
convenient. 

As in prior work on the freezing of hard-sphere mixtures [p|-p^[ , assume that the density of each species at a site is 
given by an isotropic Gaussian distribution centered at the site. (In their work on the fcc-bcc interface, Sengupta, et 
al., used a more general form of 'anisotropic Gaussian distributions' and found that there was significant anisotropy 
in the BCC density, but very little in the FCC density.) For simplicity assume that this is cut off at the Wigner-Seitz 
cell about the lattice point. (And very quickly the further assumption will be made that the Gaussian distributions 
are sharply peaked enough that the integral over the cell is equal to the integral over all space. Thus the assumption is 
essentially that the Gaussian distributions arc sharply enough peaked that there is no overlap between the distributions 
at different sites.) Referring to the cell at Ri as U then in Ui 

Pa{r) — Xai exp |— 7ai|f - -Ri| 2 j (6) 

Here, for a binary alloy, N and {A^}, {7^4 j}, and {73 i} (i — 1..N) parameterize the density ansatz, while {Xsi} is 
then fixed by the assumption of one atom per cell. For the disordered solid, X a i and j a i are independent of i, so the 
parameters are N d , X A d, jAd and 75 d- 
Define 

c a i = J drX ai exp^-j al \f-Ri\ 2 ^. (7) 

At this point the limit ^ a i — > 0, X a i —> p a a gives the uniform binary liquid. However, the states to be consider are 
those with j a i large enough that j v can be approximated by J. For large enough 7 QJ ;, 

/ \ 3 ^ 2 

Cai ~ X a i\ I . (8) 

\JaiJ 



IV. RESULTS AND DISCUSSION 



The computation of O from equation (p|) for the trial states described by the density ansatz proceeds as follows: 
Inserting the density ansatz in the first term of (||) gives: 



3/2 \ (9) 



(2) - — 

The last term depends not only on C^a(ri,i^) at the lattice points, but also near the lattice points. Assuming that 

C^(ri,f^) changes slowly enough over the length scale {l/jai) 1 ^ 2 that it can be expanded to second order, the result 
is: 



t(3 ij 

where 



(10) 



%W = cSWW.fl,'W)- (11) 
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Thus: 



+ Vf 

+e(^-CK 

+ § E E Ca j C S ij (p) 

a/3 ij 

jEE c « c ^- (z) hr + ^) f «^j(/'). ( 12 ) 



where 



C\ / / <-* L^LJ I I j I I I 



jf afiij = 2^ — ^2 — • lid) 



The equilibrium value of 7 a ^ in this approximation is: 

Upon substituting ( [l4|) into (|l^), the final term cancels with the 3/2 in the third term. This substantially simplifies 
the algebra, so minimizing with respect to the j a i this is: 



a i I \W7ai 



.3/2 



t7o 



+ ^ E E c «i c Pi c Sij(p) 

a/3 ij 

For a binary alloy c# can be eliminated. It will prove helpful to emphasize the relationship between eliminating ca 
or cb by using the more symmetrical m = ca — \- Thus m ranges from — ^ to ^. Then, for a binary alloy, 

-£{a + -) h &) + G-)-&)} 

+ iV (/J 2 + &(p)) + 7V™ al)e (£>! + 
1 W 

+ 2EE m ^4 2) w 
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where 
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Qi 

Q-2 



1 N 



i=l 



(2) 



.(2) 



'BBij(p) 



j 

\ E c aa a (p) + c ( bb a (p) + 2ci 2 l a (p) 



(17) 

(18) 
(19) 

(20) 



and additionally let: 



£>2 = 



a/3 ij 



- (D A + D B ) 



D 1 = (D A - Di 



E 



& = E C S ] 0>) = E 4 2 i « (p) + 4 2 £ a (p) - 2C^ . . ( P ) 



(21) 
(22) 

(23) 



(24) 
(25) 



Now consider a disordered state, specified by a concentration to and an overall density p. At at given temperature, 
this state will be in equilibrium for a pair of chemical potentials pa and ps- (Now buried in D±, D 2 .) For this 
disordered state, the mixed functional for the grand potential per unit volume is: 



~v~ 



1 



to In 



(7r/7Ad) 3/2 , 



1 



to In 



+ /o 



p{D 2 +g 2 (p))+pm(D 1 + g x ( P )) 
p-m 2 g 3 (p) 



,(7r/7sd) 3/2 , 



(26) 



Which can be minimized with respect to p and to to obtain the equilibrium density and concentration at a particular 
Pa and pb, or vice- versa. 

Taking the disordered state (rrid, Pd) as the reference (equilibrium) state, and minimizing equation (|26j ) to eliminate 
the chemical potentials gives: 



Di = Gi(pd) + mdGziPd) 



In [ 



m d 



{^/lAdf 12 J 

dlnjAd 



In 



+ m d ) 



di 



X^hBdf /2 

, 2 [2 



dlnjBd 



dm 



(27) 
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D 2 = Q2 (Pd) + Pd 



v 



m d ,Pd 



lad 



dp 



m d Pd 



dQx 
dp 



1 



m 2 d G?,(pd) 



2 mdPd dp 




2+ m d 

i^hAdf 12 , 
- 5(3 + "**)"* -5^- 

. 3/ X . d\n-f Ad 
+ 2^ +md ' Pd ^p— 




m d m d 



51n7 B( 



m d p d - 



dm 
dlnj Bd 



m d ,p d 



m d ,Pd 



dp 



(28) 



™d,Pd 



= /o 



2 dQ 2 
Pd ^p~ 



2 dQ x 
- mdPd -7- 

m d)Pd 



Pd 



1 2 2 



I'd 



dp 



md,Pd 



\ 2 <91 n 7sd 
™d)p d 



9p 



m d ,Pd 



6 ^ 

/3 



C/3 FapiPd), 



(29) 



Where the formula for Di has been used in computing D 2 - 

For studying the ordering transition this is the desired result for the disordered state. However, it is interesting to 
take a detour and ask what this says about the disordered state. 

One question to ask is: "How does the density vary with concentration at constant pressure?". The pressure (for 
an equilibrium state) is given by — p — Cl/V. Thus: 



Pp = fo 



2 dQ 2 
dp 



2 dQi 

m dp d ~r 

dp 



1 2 2^3 



3 (i , \ 2 d^lAd 

-liz+^Pi—gy- 

and switching back from m to c for a moment, this is: 



m d ,Pd 



-(- 

2 l2 



Tn d )p 2 d 



<9hi7_B ( 



dp 



(30) 



md,pd 



Pp = -fo 

1 



Gaj3 



+ - 2 C AP 

3 



2 dQAA 2 dQsB 

Ca ^p~ + Cb ^p~ 

2 d\nj Ad 



2cacb 



dp 



(p)- 



dpAi 
dp 



- 2 CBP 



i d In -f B < 

dp 



(31) 
(32) 



For a given p, m, and T, this equation implicitly gives the equilibrium density of the disordered solid, assuming it 
exists. 

Some insight into this equation can be gained by imagining that the Q a p are dominated by nearest-neighbor 
contributions. One then expects that for a single species the equilibrium density is give by a density near the 
minimum of Qaa- (That is, the nearest-neighbor distance in the solid is approximately equal to distance at which the 
liquid's direct correlation function has its main peak.) If the two species are significantly different in size, it plausible 
to assume that the terms involving Q a p in ([h]) are more important than those involving "f ad - Suppose the latter are 
small. At the pressure where (3p + f = 0, the density corresponding to this pressure is given by: 



G 



(33) 



And further simplifying by assuming C^ B (^0 = ( ^aa ( 





.(2), 



Cgg(r) 1/2, this reduces to 



ca — ; h (1 - c A j- 



dp 



dp 



(34) 



Considering nearest-neighbor contributions only and assuming that both C^\(r) and C BB (r) can be expanded as 
quadratics about their minima over the relevant region: 



-.(2) 



1 2 

QaA Vnn) = K A + -jk A (r nn - fA) 
1 2 

QBB{r nn ) = K B + -k B (r n „ - r B ) 
where r„„ is the nearest neighbor distance corresponding to a given density, then equation (|3^) gives: 

c A kA , (1 - CA)fc B 



CAkA + (1 - CA)fcs ^ + ca&a + (1 - CA)k 



~ r B ■ 



(35) 



(36) 



And if kA — fcs, this makes r„„ linear in concentration. 

Returning to the analysis of the ordering transition, the next step is to compute the difference in Q between an trial 
state and an equilibrium disordered state. At the transition, both the reference disordered state and the (partially) 
ordered state at the transition will be global minima of Q, and the difference in O will be zero. Hopefully the 
approximation to CI is good enough, and the range of trial states provide by the density ansatz is generous enough, 
that the solution provided will approximate the actual transition. Only trial states with 7 Q i given by the last line of 
equation ( fL6| ) need to be considered. The trial states are therefore considered as a function of p and {mi}, md, pd, 
7a d: and 7s d refer to the equilibrium state for the given chemical potentials and temperature. Combining ( |l6"| ) and 
@ through A 



A/30 « J2 \ 



2 + m i S 



IE 



- +m, 



+ m d 



In 



lAi 
7 Ad 



1 



to, 

2 



In 



In 



7b i 
7s d 



7a i 



-(N-Nd)l(l + m d )pd 



51n7A< 



m d 




m d 



d P m d , Pi 
d In 7a d 



+ (5 - m d )p . 



9 In 7b d 



dm 



m d 



dp 
dlnj Bd 



m d ,p d 



dm 



+ n (g 2 ( P ) - g 2 ( Pd )) — (n - N d ) Pd ^ 

dp 

+ Nm ave {Giip) - Q x {p d )) - (N - N d ) m dP d ^ 

dp 

1 N 

+ 2 S 51 ( TOi ~ TO<J ) ( TO J ~ TOd ) C y 2) ^ 
+ N I m a „ e 



/',/ 



m d m d (Qz{p) - Gz{pd)) 



- l - {N ~N d )mlpd d ^ 



(37) 
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The first term is the ideal mixing entropy for a pure Ising-like model. It is the same term that appears in the lattice gas 
formalism. More interesting is the second term (corrected by the next two term which cancel its linear parts) which 
also derives from the ideal part of the free energy. (*y a d/7r)~ 3 / 2 has units of volume, and represents the approximate 
volume over which the atom is likely to be. Writing this term as 



it can be thought of as representing the difference in an entropy term based on the average volume available to the 
atom to wander in. The assumption of Gaussian distributions approximates each atom as an independent oscillator 
from this point of view. We can therefor consider this expression as a lowest non-zero order approximation to the 
vibrational entropy difference. Note that this "entropy" term, which does not appear in the lattice gas formalism, 
does not directly depend on the density. Rather it depends on 7, which is the variable describing the probability 
distribution for an atom about its site, in our very simple density ansatz. This makes sense. In the lattice gas 
model there is just an occupation variable at a site, with no idea of a fluctuation of an atom position about a site. 
The non-overlapping Gaussian distributions we have adopted correspond to attaching an atom to each site with an 
independent harmonic spring. (Where the spring constant depends on both the overall density, and the state of 
chemical order.) In this picture it is not surprising that the entropy can be represented in term of the mixing entropy 
plus an entropy term based purely on the effective volume an atom occupies about its site. Thus this term represents 
a lowest-order estimate of the vibrational entropy change on ordering. Even in this simple approximation, however, 
this term depends on both the overall density and the state of chemical order, as seen below. 
The remaining terms can be rewritten as 



A/30 » . . . + ±£ ~ m d ) (m j - m d ) c\f (p) 

i=l j 

+ N (m ave - m d ) (g**(p) - G**{p d )) 

+ \n (g*{ P ) - g*{ Pd )) -{n- N d ) Pd ^- 



(39) 

Pd 



where 



g* 0) = 2g 2 (p) + 2m d g x ( P ) + m 2 d g 3 ( P ) 

= CAd GAA (p) + 2ca d C Bd gAB (p) + C 2 Bd gBB (p) 

g**(p) = gi(p) + m d g 3 (p) 

— CAd gAA (p) + {cBd~ CAd) gAB ~ CB d gBB (p) ■ (40) 

The first term is the second order term that appears in the lattice gas formulation, and the are available for 
many systems from experiment or first principles calculations. Note that it is most conveniently written in terms of 
the direct correlation function at the trial lattice constant. The last term is 0((p — p d ) 2 ) (the 0(p — p d ) portions of 
the two sub-terms cancel on expanding Q*{p) about p d ) and can be estimated from the bulk modulus of the reference 
state. The middle term is 0((m ave — m d )(p — p d )), and can be estimated from the concentration dependence of the 
density of the disordered alloy. 

Equation (|37| ) is the primary result. It is expressed in terms of the desired variables, p and {m^}. Except for the 
"vibrational entropy" terms the coefficients can be estimated using available information. It is used to predict the 
transition from the disordered to the ordered state in the same way as the lattice gas analogue of classical density 
functional theory is. Given {p d ,m d }, the disordered state always minimizes the free energy, by construction. As the 
temperature is lowered, there may be another minimum of the free energy. The highest temperature where there is 
an ordered state which minimizes the free energy, and for which the free energy difference is zero is the predicted 
transition. If the ordered state at the transition differs from the disordered state by a finite amount the transition is 
predicted to be first order. For a second order transition the ordered state will be equal to the disordered state at the 
transition, but can be determined by the fact that there will be ordered states different from the disordered states 
that minimize the free energy, and have free energies lower than the disordered state, for any temperature strictly less 
than the transition temperature. 
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One difficulty exists. The vibrational entropy term is still expressed in terms of the liquid partial direct correlation 
functions. These are unlikely to be available and suitable for many systems. It would be even more appropriate to 
use the partial correlation functions of the solid, but these are also unlikely to be available. 

What is 7 like in this theory? From equation (|lj) we have: 

k j k 



1A, 




This talks about the curvature of the liquid partial direct correlation functions. Note that 

\ dx\ ( lAp)) + nj (p)^ 



J tT L = (P)) + AtCS'to [p)) (42) 



At constant density this approximation is linear in concentration for a disordered state. The density dependence will 

(2) 

be complicated, however. Considering nearest neighbors only, we expect to be near the minimum of C a l - ,(p), subject 
to the competition between three different terms, and the contributions of the further neighbors, etc. Over a small 
range the curvature might have a simple form, but it eventually must go to zero in each direction. Because of the 
sum over the neighbor shells, the depend on the state of chemical order as well. 

In what limit does this reduce to the lattice gas expansion? For large bulk modulus the density change will be 
small. However even for fixed density, the dependence of -f a j on the ordering will appear. Thus to regain the lattice 
gas formula, an additional assumption that 7 Q i — ^ a d as well as p = p^ must be made. 



V. CONCLUSIONS 



The above computation demonstrates how a truncated binary density functional expansion has a form as the terms 
appearing in the equivalent lattice gas expansion, plus additional terms. These additional terms are the expected 
"interaction" terms from the added density variable, and an additional ideal gas entropy term involving ln(7) This 
computation is easily extended, under the same assumptions, to a global strain tensor variable in place of the density 
variable. The results are the same, except that a particular tensorial form for the "interaction" terms involving the 
traceless portion of the strain tensor is chosen from the possible forms, because the C^ 1 depend on distance only and 
have no angular part. (This will not be the case if the expansion is extended to .) 

In order to use this formalism in situations where the liquid correlation functions are unavailable some assumption 
will need to be made in regard to the l ln(7)' terms. One way to proceed is to ignore equation ( |l4| ) and make an ad-hoc 
assumption for the 7's. In a "purely harmonic" approximation we could assume 7 to be a constant, independent of 
density, concentration, and order. In this case it drops out of equation ( |37|) entirely. A slightly more ambitious 
assumption would be to assume that for the "volume" of the distribution scales as the volume per atom. In this case 
the density dependence of the ln(7) term is numerically very small compared to the other terms, so that in numerical 
solutions it would be negligible. 

This application of this formalism, including a general global elastic strain tensor rather than just the density as 
considered here, to Nickel-rich Nickel- Vanadium alloys is in progress, under the assumption that the l ln(7)' terms are 
small ||. Work is also planned to apply the formalism to the simpler case of a volume change only in CU3AU, where 
estimates of the change in 111(7) can be made from embedded atom simulations. 
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